Self-assembled granular walkers 
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Abstract 

Mechanisms of locomotion in microscopic systems are of great inter- 
est not only for technological applications, but also for the sake of 
understanding, and potentially harnessing, processes far from thermal 
equilibrium. Down-scaling is a particular challenge, and has led to a 
number of interesting concepts including thermal ratchet systems and 
asymmetric swimmers. Here we present a system which is particularly 
intriguing, as it is self- assembling and uses a robust mechanism which 
can be implemented in various settings. It consists of small spheres of 
different size which adhere to each other, and are subject to an oscil- 
lating (zero average) external force field. An inherent nonlinearity in 
the mutual force network leads to force rectification and hence to loco- 
motion. We present a model that accounts for the observed behaviour 
and demonstrates the wide applicability and potential scalability of 
the concept. 

The fascination with small-scale locomotion systems in biology, such as 
ciliates, flagellates, and molecular motors, has spurred a large number of 
attempts to construct artificial devices with similar merits [U El |3j EJ El El 
El 13 HOI ED] • Systems of particles in fluid dynamics settings can exhibit 
directed motion under periodic and symmetric forcing via the symmetry- 
breaking of surface or streaming flows fF2\ IT5] . Additionally, the controlled 
transport of particles and structures can be used for applications such as tar- 
geted delivery and stirring in lab-on-a-chip devices [T2l [TH [15] . Dry frictional 
ratchet devices composed of aspherical particles with complex interactions 
have also been shown to convert periodic and symmetric external forces into 
net locomotion [TBI IT?] . Here we report the discovery of a strikingly simple 
ratchet-like system which self-assembles in proper environments. 

When a bidisperse mixture of glass beads is moistened by a fluid and 
shaken vertically and sinusoidally, small clusters of beads occasionally take 
off from the surface of the pile and rapidly climb up the container walls 
against gravity. These clusters are held together and against the wall by 
capillary bridges; they are led by a large bead with one or more small beads 
trailing below. In this system, the self-assembly of these structures is assisted 
by the Brazil Nut Effect [18], as the large beads are transported to the top 
of the pile under vibration. Many different ascending clusters have been ob- 
served, which differed greatly in the number of involved spheres. Figure [jji 
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and Supplementary Movie Ml show a realization of such an assembly which 
spontaneously formed and climbed out of a pile of glass spheres wetted with 
a glycerol-water mixture. The upper surface of the granular pile can be seen 
as the dark region at the bottom of the images. This effect is robust, as 
it has been successfully reproduced using numerous different wetting liquids 
(silicone oil, glycerol-water mixtures, ethylene glycol), container materials 
(glass and polystyrene) and geometries (cylindrical, rectangular). 

In order to investigate this locomotion mechanism in a controlled setting, 
we have reproduced this effect with artificially assembled clusters of preci- 
sion spheres on a silicone oil wetted, chromium coated glass substrate. The 
chromium coating ensured the absence of static charging effects on the mov- 
ing clusters. The substrate was adjusted horizontally, so that the effect of 
gravity on the cluster locomotion can be neglected. A horizontal harmonic 
vibration was applied to the substrate; therefore, the horizontal acceleration 
of the substrate in the reference frame of the laboratory is a(t) = a cos(2irft), 
where ao is the peak acceleration of the substrate and / is the shaking fre- 
quency. We focussed our study on the simplest structure: an asymmetric 
dimer built out of one large and one small bead, which we call a walker. 
Quite remarkably, the walkers were found to align with the axis of vibration 
of the substrate as soon as the vibration was applied, and to migrate in the 
direction of the larger sphere. The walkers travelled with a constant speed 
when viewed stroboscopically as described in the Methods section; one ex- 
ample of this motion is shown in figure [TJd (see Supplementary Movie M2). 

Figure [2] shows our experimental measurements of the walkers' velocity 
for vibration frequencies / from 60 to 90 Hz (symbols) against the peak ac- 
celeration ao of the horizontal vibration, for a walker with R\ = 0.3 mm and 
i?2 = 0.2 mm. The values of the walkers' locomotion velocities were obtained 
by averaging a minimum of six measurements, where half of the measure- 
ments were done with the large sphere facing to the left and half where the 
large sphere was facing to the right. Both locomotion directions along the 
axis of vibration were thereby sampled, thus ruling out the possible influ- 
ence of any residual asymmetry of the setup. The size of our experimental 
errors corresponds to the difference between the largest and smallest velocity 
measurements; the experimental uncertainties represent the dispersion of the 
velocities due to the presence of dirt or surface irregularities. We observed a 
nonlinear dependence of the walkers' locomotion velocity on the peak accel- 
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eration; it rises sharply as the acceleration increases until it reaches a plateau 
near ao ~ 6 g (or 58.8 m/s 2 ). Additionally, the data exhibit no significant 
frequency dependence, which suggests that the ratio of the applied force to 
the cohesive force is a relevant control parameter (as opposed to, for exam- 
ple, an energy ratio that would imply a dependence on the peak velocity 
of the substrate [IS]). For the following data analysis, we thus introduce 
the dimensionless peak force To as the ratio between the peak inertial force 
and the average capillary force between the beads and the substrate. For a 
perfectly wetting liquid such as silicone oil, 

where M is the sum of the masses of the spheres, R is the arithmetic mean 
of their radii, and 7 is the surface tension of the wetting fluid [20]. The inset 
shows that the minimum dimensionless peak force required for locomotion 
to occur, r c , is also independent of the vibration frequency within the ex- 
perimental uncertainties. We interpret this threshold as being due to finite 
dry solid friction at the points of contact. 

We have also observed that the velocity of the walkers increases with the 
degree of asymmetry of the dimer. The asymmetry of the structure is charac- 
terised by the angle 5 between the centers of the spheres and the horizontal 
as shown schematically in the inset of figure [3j the asymmetry parameter 
is defined as sin 5 = . The velocity of the walkers as a function of 

sin 5, measured with a fixed dimensionless peak force T = 0.5, is shown in 
the main panel of figure [3j The velocity was found to vanish for symmetric 
walkers composed of spheres with identical radii (5 = 0) on a well-levelled 
substrate. We additionally observed that the symmetric dimers aligned per- 
pendicularly with the external acceleration, in contrast to the asymmetric 
walkers which self-oriented in the shaking direction. These observations es- 
tablish the asymmetry parameter sin 5 as another control parameter. 

Having identified the relevant parameters, we will now propose a simple 
model to explain our experimental observations. The purpose of our analysis 
is to show that even for the simplest assumptions to be made on the contact 
characteristics, including linearity, the asymmetry of the system inevitably 
leads to a net motion in one direction. 
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We consider a walker that is aligned with the shaking direction on a hori- 
zontal substrate. Assuming that the bead-bead and bead-substrate distances 
remain fixed, writing down the equations of motion for each bead leads to two 
sets of equations. The first one relates the magnitude of the normal contact 
forces acting at the three contact points, represented by the 3-component 
vector F := (Fx, F 2 , F 3 ), to their values at rest, represented by the vector 
F°, and to the friction forces at the three contact points F" := (f|, f|, fJ). 
Here the indices 1 and 2 refer to the contacts between spheres 1 and 2 and 
the substrate, and the index 3 refers to the contact between the spheres. We 
obtain for the normal contact forces: 

F = F° + OF", (2) 

where D is a 3x3 matrix which only depends on the angle 5 and on the rela- 
tive mass difference [i = mi ^ I m2 , where rrii is the mass of bead i [21 J . 

The second set of equations describes the rate of change of the sliding 
velocities at the contact points, Vi = Ri<pi + v x for i e {1, 2} and V3 = Rx0i + 
R2^2, where if>i is the angular velocity of bead i and v x is the translational 
velocity of the walker on the substrate: 



- GF 11 . (3) 



where G is a 3x3 matrix which depends upon the moments of inertia of the 
spheres [2T] . 

For real contacts, the friction force f} usually depends on the normal 
contact force Fj and is also an odd function of the sliding velocity Vi [221 E31 
El]: 

F}(F i ,-v i ) = -F}(F i ,v i ). (4) 

This type of friction law, associated with the relationship between the normal 
contact force F and the three friction forces described by equation [2j leads 
to a nonlinearity in the velocities in equation [3j It is this nonlinearity 
which allows for mechanical rectification effects, and thus locomotion under 
a periodic and symmetric inertial force —Ma(t). This is a generic feature 
that arises even for the simplest friction law obeying relationship |4[ a friction 





5 



force with a linear dependence on both the sliding velocity and the normal 
contact force. Indeed, if we assume a friction law f} = KV{Fi for all three 
contact points, equation [2] leads to a nonlinear expression for the friction 
forces with respect to the sliding velocities: 

F" = W(]1-DW)" 1 F° (5) 

where the 3x3 matrix W has kv\, kv 2 , and kv 3 in its diagonal as the only 
non-zero entries. 



We analytically solved equation [3] with the friction forces given by equa- 
tion [5] for a quasi-static symmetric excitation ±T /\/2, neglecting the weight 
of the beads in comparison to the capillary forces, so that at rest the normal 
contact forces F° were equal in magnitude and opposed in direction to the 
capillary forces [21J. We obtained the following expression for the average ve- 
locity of the walker v x , averaged over one symmetric excitation, as a function 
of the squared peak dimensionless force Tq and the asymmetry 8: 



where 
and 



x (Tl,6) = T ^(h 1 (Tl,6) + h 2 (Tl,6)) (6) 



cos 3 8(1 + //sin 8) 
[cos 2 8(1 + sin 2 8) - S(l + » sin 5) 2 } 2 - [2 cos 2 5 sin 5} 2 



2 4 cos 8(1 + sin 2 5) 

ri 2 {i , 0) 



[cos 6 8 - ^ (ji + sin 8) 2 } [1 + 3 sin 2 8} 



Note that if the beads were to exchange their positions, the average veloc- 
ity v x changes sign, since v x (n — 8) = —v x (6), as was observed experimentally. 

The average velocity v x (eqn. [6]) is plotted along with the experimen- 
tal data in figures [2] and [3j In order to fit this expression to our data, the 
walker asymmetry parameter sin 8 was held fixed at 0.2 in figure |2j and the 
dimensionless force T was held fixed at 0.5 in figure [3j in correspondence 
with the experimental conditions. The fit parameter k = 0.033 s/mm was 
used in both figures. While we have obtained a single well-defined value for 
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the velocity dependent friction coefficient k, we cannot compare this find- 
ing with literature values because of the unconventional friction law that 
was implemented to demonstrate the robustness of the rectification effect 
[23] . In figure [3] we observe that the predicted velocity dependence of the 
walker on the asymmetry parameter, sin 5, is consistent with the experimen- 
tally observed trend. In particular, the locomotion velocity vanishes as the 
asymmetry parameter goes to zero. In figure [2j the initial increase of the 
locomotion velocity as a function of the dimensionless force Tq is success- 
fully reproduced by the theory. However, the threshold for motion at low 
dimensionless force is not reproduced; the implementation of a more realistic 
friction law including static friction would be necessary to also predict the 
onset of motion. Another deviation is observed at higher dimensionless force, 
as the theoretical expression continues to increase while the experimental val- 
ues level off and form a plateau around Tq = 0.6. This may be qualitatively 
understood as follows. The model we proposed is only valid as long as the 
bead-bead and bead-substrate distances remain fixed. As soon as the driving 
becomes strong enough to pull the beads off the substrate or from each other, 
the model doesn't hold anymore, and another mode of motion not described 
by this model can be triggered. The plateau in the experimental velocities 
observed at T > 0.6 suggests that directed locomotion can also be obtained 
for other modes of operation than the one we described in the model. 

To investigate this in more detail, we have directly measured the rota- 
tion velocity for each sphere in a walker structure as shown in figure |4j This 
was done by tracking the angular displacements of very small tracer particles 
placed on each of the two spheres. While obtaining data using this method 
was very difficult for small dimensionless force, good quality data were suc- 
cessfully obtained for r > 0.7, for a walker composed of larger spheres. 
Figure |4^i shows that the velocity of the walker (red squares) is dominated 
by the rotation of the larger sphere 1 (black circles) at an applied acceleration 
of 2 g, T = 0.74, and a vibration frequency of 80 Hz. The smaller sphere 2 
(blue triangles) shows some sliding motion as the difference between the total 
velocity of the structure and the rotation velocity gives the sliding velocity. 
From this data we conclude that the small sphere is dragged by the larger 
sphere during each stroke of the shaker while the large sphere rotates almost 
exclusively, keeping firm contact with the substrate. In contrast, figure |4]d 
shows the walker velocity and rotation velocities of the large and small sphere 
for an applied acceleration of 10 g, Tq = 3.5, using a higher viscosity silicone 
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oil, Si AK 10 (see Supplementary Movie M3). Figure shows that there is 
considerable sliding motion present in the locomotion of both the large and 
small spheres, however the mean motion of the walker is again dominated by 
rotation of the larger sphere 1. This was determined by observing that the 
mean rotation velocity is larger (13.3 mm/s) than the mean sliding velocity 
(-5.3 mm/s). These measurements demonstrate that the walkers display dif- 
ferent modes of locomotion. 

Finally, another interesting feature is that the walkers are observed to 
remain aligned with the direction of the excitation. In order to understand 
the stability of the walker's direction of locomotion, we have to consider the 
sideways motion of the walker in the plane of the substrate. Let ip be the 
angle between the symmetry axis of the walker, projected onto the plane of 
the substrate, and the direction of the applied oscillation. If the pair of beads 
roll sideways without sliding, and the walker's acceleration associated with 
locomotion is small relative to the acceleration of the substrate, the dynamics 
of the angle ip is equivalent to the motion of a parametric pendulum. For 
small angles ip an d a harmonic driving force oscillating with a frequency /, 
the equation of motion for sideways rolling reduces to a differential equation 
of Mathieu type [25]. Defining the dimensionless time variable $ = irft, this 
equation can be written in canonical form as 

ft + [u _ AV cos(2tf)]V = (7) 

where U is proportional to the time-averaged driving force, so here U = 0. 
Our system thus dwells on the abscissa of the stability diagram in the (V, U) 
plane, whose central interval [—0.45; +0.45] lies within a stable region [25J. 
It can be shown that |V| < 0.4 for all of the parameters investigated in this 
work [21]. Since this is well within the stable region, this explains why the 
walkers remain aligned with the direction of the excitation. 

In conclusion, we have reported the discovery of a novel ratchet system 
which is notably robust, self- assembling, self-adjusting to external driving, 
functions under a wide range of experimental conditions, and displays a rich 
dynamical behaviour with several modes of motion. We have presented a ba- 
sic model of this system that successfully reproduces the main experimental 
features for moderate driving T . This model considers only the geometry 
of the system, a friction law depending both on the sliding velocity and the 
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normal contact force, and the ratio between the driving and adhesion forces, 
but doesn't rely on the exact nature of those forces nor on the average bead 
size. This suggests that similar effects could be expected for a wide variety 
of analogous systems, with different driving methods (inertial, acoustic, vis- 
cous, electromagnetic), and possibly at a much smaller scale with adhesion 
arising from van der Waals forces instead of capillary forces. 

Methods: 

The sample shown in figure [TJ left panel, consists of a wetted 50:50 mix- 
ture of glass spheres with radii ranging from 0.3-0.315 mm and 0.5-0.59 
mm, density 2.5 g/cm 3 , purchased from Whitehouse Scientific. The sam- 
ple was wetted 1% by total volume with Glycerol 87% from Merck, with 
density p = 1.23 g/cm 3 , surface tension 7 = 61 ± 1 mN/m and viscos- 
ity v = (9.8 ± 0.8) x 10 -5 m 2 /s. The sample was vertically vibrated in 
a polystyrene container with square cross-section (10.6 mm long and wide, 
44.8 mm high) purchased from VWR International GmbH. 

The remaining experiments were performed using precision ruby spheres 
with 0.2, 0.25, 0.3, 0.35, 0.4 and 0.6 mm radii and density 4.0 g/cm 3 pur- 
chased from Sandoz Fils SA. The experiments shown in figure [4] used tracer 
spheres with radii in the range of 0.038-0.045 mm purchased from White- 
house Scientific. The wetting liquids used were filtered Wacker silicone oils 
Si AK 5 or Si AK 10. Si AK 5 has a kinematic viscosity v — 5 mm 2 /s, surface 
tension 7 = 19.2 mN/m and density p = 0.92 g/cm 3 . Si AK 10 has v — 10 
mm 2 /s, p = 0.93 g/cm 3 , and 7 = 20.2 mN/m. 

Vibration was applied using a LDS model V406 electromagnetic shaker. 
We monitored the applied acceleration using a Kistler triple-axis accelerom- 
eter 8690C50 and ensured that the unwanted acceleration in the vertical 
direction was below 3% of the driving acceleration. The horizontal vibration 
experiments used either Marienfeld glass microscope slides as the substrate 
(76 mm long, 26 mm wide, 1 mm thick) or a 90 degree glass prism (with a 
length and height of 25 mm, 35.4 mm hypotenuse). A 50 mm by 25 mm area 
of the glass slides and one entire face (25 mm by 25 mm) of the glass prism 
were coated with a 4-5 nm thick layer of chromium using a BOC Edwards 
Auto 306 Evaporation System to eliminate static charge effects on the mov- 
ing clusters of beads. The substrates were wetted with a 4-6 pm thick layer 
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of Silicone oil. 

Video imaging was performed using a PCO 1200hs CMOS camera. PCO 
Camware software was used for camera control and image acquisition. Top 
view experiments (figures [2] and [3]) were imaged by mounting the camera 
onto a Zeiss Stemi 2000 C stereo microscope with a 0.63 x front lens; the 
time interval between successive pictures was set to an integer number of vi- 
brations, typically one image per shake or one image every 10 shakes, so that 
the substrate appeared immobile in the time series of pictures. Side view 
experiments used the camera with an Edmund Industrial Optics 4.5 x zoom 
lens (figure [T^i) or an Edmund Industrial Optics 10 x zoom lens (figure [1]); 
the frame rate was set to 20 images per shake. The motion of the substrate 
(figure [4]) was determined by imaging a ruler fixed to the substrate simulta- 
neously with the walker; from this the phase of the vibration was determined 
while the vibration amplitude was calculated from measurements of the peak 
applied acceleration. 

In order to obtain velocity measurements from the recorded images, the 
time series of images were opened with ImageJ software (freely available from 
|http://rsb. info .nih.gov/ij /] ) . A single threshold was applied to the images to 
isolate the structure from the background and the particle analysis package 
was used to determine the centroids of the structures in each image. Further 
analysis of the time series of structure positions was done using Matlab soft- 
ware. 
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a b 




Figure 1: Time series of ratcheting structures, a, This structure was 
self-assembled from a mixture of glass beads with radii in the size ranges of 
0.3-0.315 mm and 0.5-0.59 mm and is migrating up the sidewall of a rectangu- 
lar polystyrene container under vertical sinusoidal vibration (see supplemen- 
tary movie Ml). The sample was wetted 1% by volume with a glycerol-water 
mixture and the container was shaken at a frequency of 170 Hz and a peak 
acceleration of 16.2 g, where 1 g = 9.8 m/s 2 . The time interval between 
images is 2 s. b, This artificially assembled dimer, composed of precision 
spheres with 0.2 mm and 0.3 mm radii, is migrating along the axis of vibra- 
tion of a horizontally aligned glass microscope slide which was shaken with a 
frequency of 80 Hz and a peak acceleration of 4 g (see supplementary movie 
M2). The substrate was wetted with of silicone oil Si AK5. The time interval 
between images is 0.3 s. 
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peak acceleration a 

Figure 2: Experimental measurements of the walkers' velocity for 
varying peak substrate acceleration and dimensionless force To 
(eq. [I]). These data were obtained for a walker with Ri = 0.3 mm and 
i?2 — 0.2 mm using the following vibration frequencies: 60 Hz (diamonds), 
70 Hz (circles), 80 Hz (triangles) and 90 Hz (squares). These structures were 
assembled on a glass slide wetted with Si AK5. The solid line represents 
a fit of equation [6] to the experimental data where the coefficient of sliding 
friction, k = 0.033 s/mm, is the only fit parameter. Inset: the dimensionless 
force at the onset of motion, T c , as a function of the vibration frequency /. 
T c was measured by gradually decreasing the peak plate acceleration a and 
waiting for the walker to stop moving. The dashed line is a guide to the eye 
situated at the mean value of Tq- These experiments were performed on the 
surface of a glass prism wetted with Si AK5. The error bars represent the 
range of data resulting from multiple experimental measurements. 
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Figure 3: Quantification of the asymmetry and dependence of the 
walkers' velocity on the asymmetry. Inset: Schematic view of a walker 
composed of two spheres, 1 and 2. We characterize the walkers' asymmetry 
in terms of the angle 5 formed by connecting the centers of the spheres 
with the horizontal, and our asymmetry parameter is defined as sin 5 = 
jfripjk • The vibration axis of the substrate is denoted as x, and the angular 
displacements of the spheres at each shake are denoted as (pi and (p2- Main 
panel: The walkers' velocity as a function of the asymmetry parameter sin 5 
(squares). These measurements were performed for walkers constructed from 
precision spheres with 0.2, 0.25, 0.3, 0.35 and 0.4 mm radii on the surface 
of a glass prism wetted with Si AK5. The vibration frequency was kept 
fixed at 80 Hz and the peak acceleration was varied in order to keep the 
dimensionless force T (eq. [TJ fixed at 0.5. The error bars represent the 
range of data resulting from multiple experimental measurements. The solid 
curve shows the expected theoretical dependence of the walkers' velocity on 
the asymmetry parameter sin 6 from equation [6] with k = 0.033 s/mm (the 
same value as was used in figure [2j . 
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Figure 4: Modes of motion. The velocity (red squares) of a walker struc- 
ture in the reference frame of the substrate, where Ri = 0.6 mm and R2 = 0.3 
mm. The walker was shaken with a frequency of 80 Hz and: a, a peak sub- 
strate acceleration of 2 g, T = 0.74, using Si AK 5, b, a peak substrate 
acceleration of 10 g, T = 3.5, using Si AK 10 (see supplementary movie 
M3). The walker's velocity (red squares) is displayed alongside the rotation 
velocity —Rip of both the large (black circles) and small spheres (blue tri- 
angles) comprising the walker. The angular displacements were obtained by 
tracking the displacements of 0.04 mm radius glass spheres placed on each 
bead in the structure and were recorded at 1600 fps or 20 images/shake. 
These experiments were performed on the surface of a glass slide. The mean 
velocity of the walker is: a, 2.3 mm/s, and b, 8.0 mm/s. 
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